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Problems  in  Modeling 
the  Earth's  Trapped 
Radiation  Environment 


I.  INTRODUCTION 

The  purpose  of  this  report  is  to  highlight  for  AFGI.  the  problems  encountered 
in  modeling  the  terrestrial  trnpped  radiation  environment  for  the  past  14  years, 
and  to  indicate  steps  that  could  be  taken  to  improve  their  accuracy.  Since  this 
modeling  activity  has  been  confined  principally  to  the  energetic  electrons  and  pro- 
tons, the  following  remarks  will  apply  to  the  energy  ranges  between  40  kcV  and 
5 MeV  for  electrons,  and  between  100  keV  and  170  MeV  for  protons.  The  intent 
here  is  not  to  present  a review  of  the  characteristics  of  the  latest  models  produced 
by  the  National  Spoce  Science  Data  Center  (NSSDC),  since  these  are  amply  docu- 
mented in  the  first  12  references  listed  at  the  end  of  this  report.  However,  exam- 
ples from  these  models  will  be  utilized  along  with  other  data  to  demonstrate  the 
theme  of  this  discussion.  In  order  to  avoid  a duplication  of  the  numerous  figures 
and  tables  of  those  references  to  clarify  some  statements  in  this  report,  copies  of 
these  12  references  will  be  provided  ns  attachments.  At  the  conclusion  of  this 
report,  it  is  hoped  that  AFGl.  will  have  a better  understanding  of  the  complexity  of 
this  modeling  activity,  and  can  understand  what  must  be  done  in  order  to  signifi- 
cantly improve  existing  models  through  its  Program  to  Measure  and  Model  the 
Energetic  Pnrticle  Environment.  The  inherent  errors  in  the  models  will  be  dis- 
cussed in  general  terms  since  specific  errors  depend  strongly  on  the  particular 
uses  of  the  existing  models. 
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At  the  outset  It  must  be  recognized  that  there  are  a large  number  of  temporal 
and  spatially  varying  phenomena  which  affect  the  energetic  particle  populations 
In  the  trapping  region  of  the  earth's  magnetosphere.  Specifically,  the  past  and 
present  state  of  the  Interplanetary  medium  in  near-earth  space  (the  solar  wind 
velocity,  density,  and  temperature,  the  Interplanetary  magnetic  field  (IMF), 
including  the  7.  component,  and  passage  of  interplanetary  shocks)  play  a major 
role  in  producing  excitation  in  the  magnet ospheric  cavity  in  the  form  of  magnetic 
storms  and  substorms,  which  give  rise  to  particle  injection  and  acceleration 
mainly  through  the  plasma  sheet  in  the  geomagnetic  tali.  Boundary  currents  in 
the  magnetopause  produce  alterations  in  magnetic  field  within  the  cavity  ns  well  ns 
do  the  injection  of  lower  energy  plasma  from  the  ionosphere,  the  geomagnetic  tall, 
and  the  interplanetary  medium  to  form  the  ring  current  in  the  spatial  region  between 
2 and  f»  earth  rndii.  In  addition,  electric  fields  produced  by  convection  of  IMF  past 
the  earth  nnd  through  plasma  instabilities  within  the  cavity  affect  the  motion  of 
lower  energy  particles  and  contribute  to  acceleration  anti  diffusion  of  the  more 
energetic  particles.  The  literature  is  rich  in  the  description  of  these  many  com- 
plex phenomena  nnd  even  a brief  review  of  them  is  beyond  the  scope  of  this  sum- 
mary. A copy*  is  also  attached  since  it  Is  a recent  survey  of  the  trapped  radiation 
by  the  modeling  group  at  NSSIK'  and  contains  references  to  the  recent  literature  on 
this  subject. 

We  will  confine  our  discussion  to  the  following  quantities.  The  basic  measure- 
ments obtained  by  energetic  electron  and  proton  instruments  flown  on  satellites  are 
either  the  omnidirectional  or  unidirectional  flux  (pitch-angle  distribution)  within 
certain  energy  intervals  or  above  nominal  threshold  energies.  In  the  strict  sense, 
the  omnidirectional  flux  is  the  number  of  particles,  of  specified  type  and  energy 
range,  that  pass  through  a spherical  surface  with  a unit  cross-sectional  nren  in  a 
unit  time  nnd  Is  expressed  most  typically  ns  the  number  of  electrons  (or  protons)/ 
cm  -sec  wl‘h  the  appropriate  qualifiers  on  the  energy  range  observed.  Most 
practical  detectors  fall  short  of  this  Ideal  geometry  but  differ  from  the  more 
directional  instruments  which  nre  sensitive  to  particles  only  in  a narrow  angular 
range.  For  purposes  of  this  report  we  consider  directional  detector  systems  as 
those  with  a solid  angle  of  less  than  0,  0»>  sterndinns,  while  omnidirectional  detec- 
tors ar»  those  sensitive  over  substantial  fractions  of  the  4»  sterndinns  of  the  sphere. 
Because  a spacecraft  nnd  its  systems  provide  a significant  amount  of  shielding,  the 
most  ideal  geometry  for  such  omnidirectional  measurements  are  2»  strrndtan 
detectors.  Because  of  the  nature  of  the  trapped  radiation,  this  geometry  is  equiva- 
lent to  the  spherical  4»  steradian  detectors.  Directional  trapped  radlntion  detectors 

1.  Chan,  K.  W. , Sawyer.  D.M.,  and  Vette,  J.I.  (1977)  The  trapped  radiation 
population,  in  The  Trapped  Hndtntion  Handbook,  Kds.  J.B.  Cladis,  G.  T. 
Davidson,  nndT.  L.  Newkirk  DNA  2524  Change  5,  Defense  Nuclear  Agency. 
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must  be  referenced  to  .some  spatial  axis  which  is  usually  the  actual  magnetic  field 
observed  by  a magnetometer  on  the  spacecraft.  The  spatial  position  of  either  type 
of  instrument  must  be  known  during  the  accumulation  of  particle  counts  over  a 
small  period  of  time,  and  tins  requires  the  use  of  some  reference  coordinate  sys- 
tem. Consequently,  obvious  independent  variables  of  the  flux  function  are  the 
positional  coordinates,  the  energy  coordinate,  an  angular  coordinate  for  direc- 
tional systems,  and  time,  since  it  is  known  from  observation*  that  time  variations 
are  very  im|«ortant.  It  is  the  variation  of  the  observed  fluxes  with  these  variables, 
coupled  with  the  physical  characteristics  of  the  orbiting  instruments  (particle 
energy  response,  angular  response,  and  dvnamic  range  for  counting),  that  will  be 
used  to  describe  the  inherent  problems  in  obtaining  suitable  models  of  the  ener- 
getic trapped  radiation  environment.  The  extremely  large  volume  of  space  which 
must  be  observed  is  the  dominant  factor  responsible  for  the  poor  sampling  that  has 
characterized  the  basic  data  used  to  produce  existing  models. 


1 UMWIUNATI  SYSTIMS  \\U  SPATIAL  t.KAIHINTS  Ol 
r YKTia.fr  I’M  \fr> 

Although  the  trapping  region  of  the  magnetosphere  is  a time  varying  volume  of 
space,  reasonable  time -averaged  limits  of  the  region  are  within  1 1 It  of  the  earth 
and  restricted  to  t?0°  in  geomagnetic  latitude.  Since  the  energetic  trapped  particle 
represent  a small  energy  density  relative  to  the  energy  density  of  the  magnetic 
fields  within  this  region,  the  particle  motions  are  confined  to  gyration,  Knitter,  and 
drift  on  these  magnetic  field  lines.  In  addition,  the  density  of  the  earth's  atmos- 
phere restricts  energetic  particle  populations  to  the  region  above  1.  la  It  ; the 
radial  gradient  become*  very  steep  tv  low  l.  2 ll  . 

To  appreciate  the  vastness  of  the  trapping  region,  consider  that  this  volume 

must  be  broken  into  cubical  cells  of  0.  I It  on  a side  at  radial  distances  between 

e 

l .2  and  7 It  , and  0.2  It  between  7 and  1 1 It  , in  order  to  properlv  account  for 
r c o • • * 

the  spatial  gradients  of  the  observed  time-averaged  fluxes.  The  reason  for  such 

cell  sizes  will  become  apparent  later  in  the  report.  A simple  calculation  reveals 

that  if  a strictly  geographic  coordinate  system  is  used,  1.0  \ I0h  such  cells  are 

Involved.  To  provide  a reasonable  energy  spectral  distribution  about  12  electron 

and  20  proton  energies,  must  be  used.  Since  the  observed  omnidirectional  fluxes 

•» 

are  observed  to  vary  from  cosmic  rav  levels  of  order  l cm  -sec  to  order 
8 2 

10  cm  -sec  for  the  lowest  energy  particles  considered,  a straight  forward  geo- 

n l| 

metrical  approach  might  require  a storage  sire  of  some  10*  to  10°  bytes.  This 
would  clearly  be  cumbersome,  if  not  totally  impractical,  even  with  present  day 
computer  technology. 


Consequently,  the  modelers  must  use  some  physical  understanding  of  the 
trapping  phenomena  as  well  ns  data  compression  techniques  in  order  to  produce  r 
model  that  can  bo  used  in  a practical  manner.  The  first  trapping  region  of  the 
magnetosphere  that  was  understood  was  the  inner  zone  or  Van  Allen  belt.  Kxcept 
for  perturbations  caused  by  extremely  large  magnetic  storms  or  bv  artificial 
injections  from  nuclear  detonations,  the  ordering  of  particle  fluxes  in  this  region 
was  shown  to  be  adequate  by  using  the  two-dimensional  11  (scalar  magnetic  intensity) 
and  l parameter  system  introduced  by  Mcllwain  and  computed  using  the  harmonic 
representation  of  the  internal  geomagnetic  field.  Other  coordinate  variables  de- 
rived from  these  two  have  also  been  useful  and  in  the  practical  models  produced, 
the  ratio  of  the  magnetic  field  at  the  point  of  measurement  to  that  at  the  geomag- 
netic equator  of  the  same  field  line,  II  'It  along  with  l,,  have  been  used  as  the 
basic  system.  The  great  advantage  of  such  a system  is  that  it  is  two-dimensional 

instead  of  three-dimensional;  this  alone  results  in  a reduction  of  the  number  of 

2 

spatial  cells  required  by  more  than  a factor  of  10“.  In  addition,  the  use  of  It / It  ^ 
provides  a rapid  recognition  of  positions  near  the  geomagnetic  equator  on  each  l 
shell. 

In  the  gross  sense  the  major  error  produced  by  the  choice  of  a two-dimension- 
al system  is  at  very  low  altitudes  for  electrons,  where  longitudinal  asymmetries 
arr  observed  because  of  the  displacement  of  the  dipole  moment  from  the  center  of 
the  earth,  and  of  the  electron  precipitation  produced  by  atmospheric  scattering. 
However,  flux  levels  at  altitudes  below  400  km  are  small  enough  that  only  highly 
radiation  sensitive  devices  such  as  man  and  film  require  additional  concern  and. 
even  here,  low  altitude  orbit  satellites  perform  an  average  over  longitude  m 
several  days. 

As  one  moves  to  the  outer  electron  zone,  it  in  apparent  that  a third  spatial 
coordinate  is  needed  for  the  region  above  5 even  for  long-term,  time-averaged 
fluxes.  Although  it  is  understood  distortions  of  the  internal  geomagnetic  field  are 
produced  by  the  ring,  magnetopause,  and  geomagnetic  tail  currents,  a more 
physically  meaningful  coordinate  system  using  these  external  fields  is  not  practical 
nor  were  there  good  external  field  models  available  m the  early  days,  of  modeling 
this  environment.  In  addition,  nonadiabatic  time  variations  of  the  particle  fluxes 
(discussed  later)  are  known  to  dominate  the  region  above  2 II  . Continual  and  dis- 
crete flux  increases  following  magnetic  storms  to  distances  ns  low  as  1.5  |{ 
p eclude  the  use  of  external  field  models  to  order  the  trapped  fluxes  effectively. 

The  third  spatial  coordinate  adopted  for  model  work  has  been  local  time,  since  it 
Is  the  best  one  to  account  for  the  distorted  magnetic  field  with  respect  to  the  sun- 
ward and  anti-sunward  asymmetry,  which  alters  particle  motion  over  a much  larger 
volume  of  space  than  that  where  longitude  is  n better  third  spatial  coordinate.  Since 
the  dominant  use  of  the  trapped  radiation  models  has  been  for  satellites  which  orbit 


for  long  enough  periods  to  effectively  hvernge  both  time  and  local  time  variations, 
all  practical  models  have  contained  local  time-averaged  fluxes.  This  results  in  a 
two-dimensional  spatial  matrix  for  storing  the  models  in  a digital  computer  for  all 
parts  of  the  trapping  region. 

To  demonstrate  some  of  the  problems  related  to  the  spatial  flux  variations  and 
the  coverage  of  the  trapping  region  by  various  satellites,  we  employ  H /Ho.  L,  and 
the  local  time,  ♦.  as  the  spatial  coordinates.  Some  of  the  graphs  referred  to  in 
the  references  use  11  instead  of  11/11  but  this  should  not  present  any  problem  to 
the  readers.  The  omnidirectional  time-averaged  flux  is  then  a function  of  these 
spatial  coordinates  and  of  the  energy  of  the  particle  namely,  J(ll/B  , L,  ♦.  K). 
llecause  of  the  variations  of  the  flux  with  energy,  it  is  necessary  to  consider  a 
number  of  energies  (-»  10)  to  provide  numerical  models  which  can  then  be  inter- 
polated adequately  to  obtain  the  energies  of  interest  for  a given  application.  The 

flux  on  a given  field  line  (l.  shell)  varies  from  some  peak  value  nt  11  11  1,  the 

8 3 ^ 

geomagnetic  equator,  in  the  range  of  10  to  10  depending  on  particle  and  energy, 

2 

down  to  about  1 particle  ’em  -sec  at  the  atmospheric  cutoff.  With  such  large 
fluxes,  it  is  more  appropriate  to  express  radial  gradient  in  the  form  of  the  log- 
arithmic gradient,  for  example 


<L 

ill 


(In  .!> 


I dJL 

J dl. 


For  inner  rone  electron  fluxes  above  1 MeV  near  the  inner  radial  boundary, 

this  logarithmic  gradient  is  typically  ll'0/R  . Consequently,  choosing  a cell  sire 

with  a radial  dimension  of  0.  1 can  still  produce  large  errors  in  computing  the 

omnidirectional  flux  from  a numerical  model  with  fluxes  given  at  various  matrix 

points.  At  1.  values  around  2.0,  this  gradient  is  typically  5/Ht>.  To  obtain  a 

reasonable  grid  sire  for  the  models,  in  view  of  radial  gradients,  tlie  1,  space  is 

tvpicallv  divided  in  0.  I It  segments  from  I,  1.2  to  7.  0 and  to  0.  2 It  for  l.  values 
• e e 

alx've  7.0.  The  11  11  intervals  are  chosen  ' > considering  the  omnidirectional  flux 

o • * 

at  the  equator  and  assigning  aKnit  three  11/11  Intervals  per  decade  of  flux.  Con- 

7 o 

sequent ly,  if  the  equatorial  flux  is  10  particles  cm  -sec  for  a given  1.  shell,  about 
21  intervals  w ill  be  chosen.  I sing  this  approach,  the  total  number  of  11/11  , I.,  and 
K grid  points  in  the  models  is  about  7,  000  for  both  particles.  Assuming  4 bytes 
are  necessary  for  the  flux  at  the  grid  (Mints  the  models  are  stored  in  a Knit 
28,  000  bytes. 

With  the  exception  of  the  region  covered  by  geostationary  equatorial  satellites, 
it  is  very  difficult  to  obtain  a measure  of  the  local  time  variation  of  the  trapped 
flux.  The  reason  for  this  difficulty  is  that  the  sampling  of  all  values  of  local  time 
require  several  months  to  one  year  in  most  orbits.  Consequently,  the  removal  of 


Inherent  time  variation*  through  averaging  is  not  easy  and  local  time  variation* 
can  still  be  masked  by  these  time  variations.  However,  a geostationary  space- 
craft samples  all  values  of  local  time  within  34  hr  and  so  the  local  time  variations 
can  be  ascertained  fairly  accurately.  The  results  obtained  by  t'aulika*  and  ltlake 
on  ATS  l are  shown  in  Figure  33  (p  HO)  for  various  energy  electrons.  The  ratio 
of  maximum  to  minimum  varies  from  about  l.  7 at  300  keV  to  3.  3 at  1.9  MeV  for 
this  data.  The  results  obtained  in  the  construction  of  outer  rone  electron  model 
AF-4  Indicate  that  these  diurnal  variations  are  not  extractable  below  l.  • 5.0  and 
tend  to  build  up  at  least  until  one  get*  to  t,  • #.  5.  It  is  difficult  to  determine  this 
variat  on  at  higher  I.  values  and  the  average  variation  can  get  as  large  as  a factor 
of  15. 

I'ntll  more  definitive  results  can  be  determined  by  much  greater  coverage, 
the  averaging  over  local  time  is  the  only  practical  approach  for  presenting  a useftil 
model.  Most  satellites  sample  all  values  of  local  time  during  live  course  of  orbit- 
ing for  several  years,  consequently  a local  time  averaged  flux  is  the  appropriate 
form  to  be  applied  for  accumulated  fluence  determ Inatlon.  However,  in  comparing 
instantaneous  fluxes  from  two  measurements  that  differ  only  in  their  local  time 
variable,  large  differences  can  be  expected. 


3.  SATKUJTK  COVKMUK  OF  THK  TO  MT1M'.  MUION 

The  volume  of  the  trapping  region  is  so  large  that  adequate  coverage  to  moni- 
tor energetic  particle  populations  has  never  been  practical.  Since  each  satellite 
is  confined  to  a particular  orbit,  no  one  satellite  has  ever  provided  a good  sampling 
of  all  the  regions  for  which  mersurements  should  be  made.  Hecause  of  the  general 
character  of  the  t rapped  radiation,  the  most  complete  coverage  by  a single  satellite 
would  be  roe  in  an  equatorial  orbit  with  an  apogee  and  perigee  that  covered  the 
trapping  region;  such  a spacecraft  would  have  to  have  a unidirectional  particle 
detection  system  that  was  capable  of  measuring  the  complete  equatorial  pitch  angle 
distribution  over  the  desired  range  of  energies, 

To  demonstrate  that  the  lack  of  spatial  coverage  for  any  specific  epoch  is  a 
very  real  and  practical  problem,  a series  of  figures  have  been  prepared  to  show 
this  coverage  for  the  current  model*  that  have  been  produced  by  NSSIX'.  These 
current  model*  are  summarised  in  Table  1,  Proper  documentation  for  AKl-7  has 
not  been  issued  although  the  computer  decks  have  been  distributed.  A discussion 
of  the  need  for  this  latter  model  on  an  Interim  basis  (hence  the  use  of  the  I in  the 
model  designation!  will  be  discussed  in  Section  S, 

2.  Vette,  J.l. , Kd.  (19711  Model*  of  the  Trapped  Kadiation  Knvlronmont- 
Volume  Vllt  Long  Term  Time  Variations^  'KAlxA  'TT^T'ITi 
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Tafcl*  1.  Currrni  N SvX.  KadiaHon  Knvironrner.t  Models 


Figures  1 and  2 show  tho  H/ll  ' H coverage  of  the  data  utilized  in  construct- 
ing tlio  electron  an«l  the  proton  models,  respectively.  The  magnetic  eipiator  is 
represented  by  the  axis  of  It/ll  1 or  X 0°.  The  entire  area  ineluded  by  the 
100  kin -cutoff  curve  and  the  eipmtorial  axis  Is  accessible  to  spaceltorno  oxpert- 
inents.  Tho  regions,  in  which  datn  from  only  one  spacecraft  and  froi.'  ..lore  than 
one  spacecraft  were  available,  have  boon  shaded  separatelv.  The  blank  area  tndi- 
cates  no  data  coverage,  and  in  'he  case  of  A i’ll  (Figure  2)  one  should  note  th;-' 

I,  (1.  li  is  the  model  cutoff. 

With  some  2-1  spacecraft  used  for  Al’ll  and  l.'t  for  AK-J/-5  models,  and  con- 
sidering tiie  accumulation  of  coverage  over  a |>oriod  of  10  years  for  all  energies, 
it  is  surprising  to  note  that  there  are  still  gaps  ami  regions  not  well  monitored.  If 
we  had  displayed  the  cross-section  of  a month  of  time  and  a given  energy  interval, 
the  spatial  coverage  would  have  been  far  from  complete. 


Figure  I.  H/lt  - I.  Coverage  for  Kxperiinents  Itsed  in  Constructing  the  AH- 4 
anil  AK-ft  Models.  The  area  shaded  for  one  satellite  is  covered  by  the  j*olar 
orhiter  Injun  :i 
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Kllturr  12.  II  H(l  - I fovri'nisc  for  I'xprrlmrntn  t’nrd  hi  Count  rurtlnjj  llir  A 1 * -H 
Modrl.  Tlir  modrl  and  tlir  data  iwed  have  n nitoff  at  I t«.(* 


A nummary  of  tin'  rnrruy-thnr  rovrraijr  of  itata  u.nrd  In  tlir  varloim  tnodrln  In 
altnwn  In  Kl(5iirrn  1 mul  0.  In  oarli  of  tbrnr  diagram*,  tlir  rnrrjsy  intrrvala  of  tlir 
parttrlr  rlinmtrla  arc  markrd  by  tlir  vrrttral  k-irn,  ami  tlir  nominal  energy  <'f  tlir 
thrraliold  rliannrl  la  Indlcatrd  by  a nln^lr  upward  arrow.  Tlir  thnr  apan  i*f  tlir  data 
rliannrl  nard  In  iiirdrling  In  rrprrnrntrd  l>v  tlir  Irnittli  of  tlir  Imrlrontal  bar  tojjrtlirr 
with  tlir  apaorrraft  dmlj’iiallon.  Itv  Inaprrtlon  of  llir.nr  ft^irr.n.  It  In  liiiiurdlatrlv 
apparrnt  tliat  tlir  hl|{lt*rnrri:y  data,  '•'J  MrV  for  rlrrtronn  and  '100  MrV  for  proton*, 
arr  lartfrly  larking  hi  tlir  data  l*anr  of  rurrrnt  tnodrln,  I'urtlirrmorr,  tlirrr  havr 
brrn  vrrv  llttlr  data  avallal'lr  brtwrrn  10, Ml  and  101*1.  I'l^urr  t>  In  tlir  plannrd  data 
rovrraur  to  l>r  unrd  for  Modrl  A II *7  ami  will  l»r  dl.nru.nnrd  In  drtall  In  Srrtlon  (I, 
Itrturnlng  to  tlir  In-nltu  niranurrmrntn  by  an  Individual  aatrllitr,  tlir  namplltij; 
limit  of  a nl vrn  orbit  ran  br  rradllv  illuntratrd,  For  tlir  |xdar  orbitrra  with  aivgrr 
atvut  4,000  km  (for  rxamplr,  OV1-1,  A /.l  HI,  tlir  r>|uatorlat  rrjjb'ti  of  all  I ' 2,0  i.n 
not  rovrrrd.  An  rxamplr  of  tlir  artual  data  nnpiirrd  by  A /.I’ll  at  tlirrr  I -valura 
liavr  brrn  aluiwn  In  Flgurr  41  (p  7.'l),^  It  In  obvioun  that  tlirnr  |*olar  orbitrra 


1.1 


ItM  ItW  IM1  IM)  (W)  IM4  \tm  )«M  tw  IN*  t«M  w* 

VIA* 


Figure  3.  Knergy  - Time  Coverage  for  Kxperlmen's  Used 
in  Constructing  the  AK-4  Model.  The  nominnl  energy  of 
threshold  detectors  are  shown  with  the  upward  po.nting  arrows 
and  the  energy  bands  of  the  interval  detectors  are  pi-esented 
by  the  vertical  bars.  The  time  span  of  the  data  used  is  in- 
dicated by  the  length  of  horizontal  bars.  The  Starfish  Deto- 
nation event  is  marked  on  the  time  nxts 


can  contribute  only  about  half  of  the  11/11  range  required  by  the  model  and  no  equa- 
torial measurement  can  bo  made  for  L values  above  2.0. 

On  the  other  hand,  the  high-alt itude,  low-to-mid-latitude  satelltt?s,  such  ns 
000-3  and  000-5,  are  restricted  by  the  poor  time-coverage  of  the  trapp'd  radia- 
tion zones.  With  an  initial  apogee  of  about  23  It  these  spacecraft  can  sample  a 
given  L-value  only  twice  every  two  days.  This  kind  of  data  sampling  will  miss 
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twi  1MJ  «WJ  IW  W 1IM  Ml  I Mi 


VI  AN 

Figure  4.  Energy  - Time  Coverage  for  Experi- 
ments Used  In  Constructing  the  AK-ft/-fi  Models. 

Notations  are  the  same  as  In  Figure  3 

most  of  the  predominant  short-term  temporal  variations  In  the  outer  zone,  which 

will  lie  further  explored  In  the  following  sections, 

Detailed  spatial  and  spectral  covernge  have  been  given  In  most  of  the  model 

3 

documents.  Headers  can  refer  to  Figures  Ofi  through  128  (pp  U8-I4t>)  for  specific 
examples.  General  comments  on  the  near-term  outlook  on  this  area  and  recom- 
mendation for  future  measurements  will  be  dtacussed  In  Sections  fi  nnd  7. 


3.  Sawyer,  D.  M.,  nnd  Vette,  J.l.  (lP7ti)  A l’-H  Trapped  Proton  Knvlronment  for 
Solar  Maximum  nnd  Solar  Minimum.  NSSDC  '\VDC-A-lt&S  711-01'^ 

1ft 


Figure  6.  Energy  - Time  Coverage  for  the  New  Data  to  be 
Used  in  Constructing  the  AE-7  Model.  Notations  are  the  snme 
as  in  Figure  3,  Since  the  AE-7  model  has  not  been  completed, 
this  figure  Is  only  a tentative  graphic  illustration  of  a plnn 


4.  SPECIAL  VARIATIONS,  THE  CALIBRATION  OF  DETECTORS,  ANI) 

OTHER  DETECTOR  PROBLEMS  IN  MODELING 

The  large  spatial  gradients  encountered  throughout  the  trapping  region  make 
the  coverage  problem  a difficult  one  since  an  impractlcally  large  number  of  orbits 
would  normally  be  required.  Spectral  variations  produce  another  problem  that 
arfects  the  modeling  efforts.  The  majority  of  the  data  used  in  producing  trapped 
radiation  models  has  been  obtained  from  scientific  satellites  in  which  the  emphasis 
has  been  on  understanding  the  dominant  physical  phenomena  in  any  given  region  at 
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a given  time.  Surprisingly,  the  quantitative  response  of  the  instrument  to  the 
trapped  particles  has  not  been  the  important  factor  in  understanding  the  physics  of 
region;  so  whether  an  electron  detector  is  predominately  sensitive  to  all  electrons 
about  1.  G or  1.9  McV  makes  little  difference  in  the  ensuing  data  analysis  and 
physical  interpretation.  However,  the  average  energy  spectrum  in  the  equatorial 
geostationary  region  gives  approximately  a factor  of  4 difference  in  the  omnidirec- 
tional flux  for  these  two  energy  levels.  In  using  data  from  many  satellite-borne 
experiments,  a quantitative  comparison  of  different  detector  systems  is  an  impor- 
tant factor  in  obtaining  the  best  value  of  the  flux.  This  requires  that  the  energy 
sensitivity  of  a detector  system  be  known  with  reasonable  accuracy.  Since  elec- 
trons scatter  quite  readily  in  the  shielding  material  employed  with  these  detectors 
ns  well  as  with  spacecraft  mass,  the  energy  response  can  best  be  obtained  through 
extensive  calibration  employing  sources  of  known  intensities  and  energies.  There 
have  been  many  electron  detector  systems  flown  which  have  never  been  properly 
calibrated.  This  results  in  considerable  error  since  the  spectral  distribution  of 
trapped  electrons  varies  greatly  in  time  and  space.  Electron  detectors  that  meas- 
ure energies  mainly  below  1 McV  can  be  calibrated  by  many  available  sources  nnd 
accelerators,  nnd  even  reasonable  estimates  can  be  made  bv  calculations  if  the 
electronic  discriminators  are  known  to  be  set  at  reasonable  levels  to  detect  most 
of  the  energy  lost  by  the  electron  in  the  sensitive  volume  of  the  detector.  However, 
for  energies  above  1 MeV,  a detailed  calibration  is  nn  absolute  necessity.  The 

4 

ATS  6 system  of  Paulikas,  Hlake,  and  Imamoto  referred  to  in  Figure  7 certainly 
meets  the  level  of  detail  required  and  covers  an  adequate  energy  range.  We  have 
encountered  a recent  example  where  comparisons  with  these  ATS  f>  detectors  nnd 
others  claiming  nn  effective  threshold  of  2 MeV  show  a factor  15  discrepancy  in 
the  average  flux.  It  has  been  our  belief  that  most  energetic  electron  detectors 
above  1 MeV  have  never  been  adequately  calibrated;  therefore,  there  is  still  con- 
siderable controversy  about  the  best  flux  values  in  this  energy  region.  Proton 
detector  systems  have  usually  mensured  particles  within  an  energy  Interval  that  is 
too  large  for  it  to  be  considered  a differential  energy  analyzer.  Although  many  of 
these  systems  have  not  been  properly  calibrated  in  the  past,  a calculated  energy 
response  is  generally  more  accurate  than  for  electrons.  Consequently,  the  com- 
parison of  various  instruments  has  been  somewhat  more  favorable  although 
extensive  calibrations  arc  still  necessary  to  determine  flux  levels  to  accuracies 
better  than  a factor  of  two. 


4.  Paulikas,  G.A.,  Blake,  J.B.,  and  Imamoto,  S.  S.  (1975)  ATS-6  Energetic 
particle  radiation  measurement  of  synchronous  altitude,  IEEE  Trans,  on 
Aerospace  and  Electronic  Systems  AES-11 :1 138-1 144. 
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Figure  7.  Energy  Calibration  Data  for  the  ATS  fi  Detectors.  K2.  E3,  and  E4 
are  threshold  detectors.  El  is  interval  detector  (After  Paulikas  et  al.  Ref.  14) 


Returning  to  electron  detectors,  it  has  been  our  experience  that  the  interpre- 
tation of  threshold  detectors  by  various  investigators  in  terms  of  quantitative  flux 
levels  is  somewhat  inconsistent.  This  inconsistency  can  range  up  to  a factor  of 
four  in  terms  of  the  absolute  value  of  flux  in  certain  regions  of  space.  However, 
whenever  the  energy  response  curve  has  been  available  to  us,  this  difficulty  has 
been  removed  by  folding  the  final  model  spectrum  with  the  response  curves  in  a 
consistent  manner  to  determine  what  the  proper  conversion  factor  from  counting 
rate  to  flux  should  be.  For  the  higher  energy  threshold  detectors  this  factor  is 
spectrum  dependent,  therefore,  it  is  also  spatially  dependent.  Of  course  we  have 
had  no  valid  way  of  establishing  the  accuracy  of  the  energy  response  curves  pro- 
vided along  with  the  particle  data.  Therefore,  the  final  flux  values  for  the  models 
are  a subjective  best  estimate.  However,  it  has  been  our  practice  to  illustrate 
fairly  completely  the  differences  between  the  experimental  flux  values  obtained  by 
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each  experiment  and  the  model  values.  A scanning  of  the  most  of  the  references 
cited  will  demonstrate  this  detailed  comparison. 

Directional  detectors  have  presented  several  difficulties  for  incorporating  the 
results  into  the  models.  These  instruments  have  either  provided  the  flux  perpen- 
dicular to  the  magnetic  field  line,  or  given  a measurement  at  a particular  pitch 
angle  at  some  point  on  the  field  line.  In  order  to  obtain  the  omnidirectional  flux, 
the  complete  equatorial  pitch  angle  distribution  on  each  field  line  must  be  deter- 
mined. Consequently,  each  pitch  angle  measurement  must  be  translated  to  the 
equator  using  adiabatic  theory.  A typical  equatorial  pitch  angle  distribution  is 
shown  in  Figure  8.  Near  the  atmospheric  cutoff  one  can  see  that  the  slope  be- 
comes quite  steep.  This  results  in  various  errors  because  the  finite  angular  re- 
sponse of  the  instrument  and  the  direction  of  instrument  relative  to  the  field  line 
must  be  known  fairly  well.  During  injection  events  when  there  are  large  increases 
in  the  flux,  this  pitch  angle  distribution  is  known  to  change  from  the  average  or 
equilibrium  values.  Since  one  is  effectively  performing  a different  type  of  time 
averaging  than  with  omnidirectional  detectors,  the  final  values  might  be  different 
by  as  much  ns  a factor  of  two  And,  of  course,  it  is  not  possible  to  obtain  an 


Figure  8.  Kquatorial  Pitch  Angle  Distribution  of  822  keV  Klectron  Data  From 
OGO-5.  The  polynomial  curve  is  a best  fit  to  the  inner  rone  data  base  at  l.  1.  5. 
The  slopes  of  the  gradients,  l/j  Aj  An,  have  the  unit  (deg**) 
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instantaneous  comparison  of  an  omnidirectional  detector  and  a directional  detector 
because  they  are  measuring  different  portions  of  the  particle  population.  Direc- 
t tonal  detectors  require  that  the  angular  calibration  be  done  more  accurately  than 
for  an  omnidirectional  detector,  and  it  has  been  our  experience  that  the  geometric 
factors  of  such  systems  as  quoted  by  experimenter  can  change  by  a significant 
factor  even  after  we  have  been  processing  the  data  for  lengthy  periods  in  the  pro- 
duction of  models.  Generally,  the  geometric  factors  of  omnidirectional  systems 
can  be  calculated  to  better  than  a factor  of  two,  so  the  nu>deler  has  some  check  on 
the  quoted  values,  where  as,  the  calculation  of  the  effective  geometrical  factor  for 
directional  systems  require  a much  more  thorough  knowledge  of  the  construction 
of  the  instrument  and  do  depend  on  the  actual  pitch  angle  distribution  being  observed. 
Directional  detectors  are  extremely  useful  for  scientific  purposes  but  for  the  con- 
struction of  time-averaged  omnidirectional  models,  omnidirectional  detectors 
present  fewer  pivblems.  In  addition,  the  effect  of  backgrounds  produced  by  ener- 
getic particles  ttud  penetrate  the  shielding,  which  define,  in  many  cases,  the 
sensitive  aperture,  are  not  always  able  to  In*  understood  and  removed  properly. 

Other  problems  that  one  must  be  aware  of  in  using  the  data  from  particle 
detectors  are  the  saturation  effects.  Since  these  effects  can  bo  produced  by  dead 
time  within  the  sensor  device  such  as  in  Geiger  tubes,  or  in  the  discriminators 
that  determine  the  minimum  electrical  pulse  that  w ill  be  counted,  or  in  the  counters 
or  amplifiers,  it  is  difficult  to  understand  the  various  ways  in  which  saturation 
effects  manifest  themselves.  The  valid  dynamic  range  of  an  instrument  system 
are  usually  not  clearly  defined  to  the  modeler,  and  must  be  inferred  by  studying 
the  final  flux  values  or  pitch  angle  distributions  involved.  As  part  of  the  calibration 
procedure  of  an  instrument,  the  system  should  be  exposed  to  large  enough  intensi- 
ties of  various  energy  particles  to  drive  it  into  saturation.  The  saturation  charac- 
teristics should  bo  made  available  as  part  of  the  supporting  documentation  so  that 
the  modeler  is  not  trying  to  incorporate  measurements  which  are  in  gross  error 
because  of  these  effects.  We  can  cite  some  very  subtle  effects  that  we  have  ob- 
served in  data  where  the  experimenter  maintains  the  Instrument  was  performing 
properly.  It  is  very  unwise  to  base  a model  on  the  data  from  one  instrument:  but 
because  of  the  incomplete  coverage  in  the  four -dimensional  space,  H/H  I.,  K, 
time,  we  have  had  to  do  this  in  certain  regions. 


X.  TKMPORAI.  VARIATIONS  OF  Till'.  TRAITFII  T ARTM^U-.  POHLATIONS 

The  most  difficult  problem  encountered  in  producing  models  of  the  energetic 
particle  environment  results  from  the  wide  variety  of  time  variations  that  are 
observed  throughout  the  radiation  belts.  We  wll’  discuss  these  various  changes 
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briefly  because  it  is  essential  to  understand  why  they  pose  such  a problem,  and  it 
is  useful  to  know  the  approaches  that  have  been  used  in  our  modeling  efforts.  Since 
the  most  prevalent  use  for  these  models  has  been  to  determine  the  approximrte 
radiation  environment  that  new  missions  will  encounter  in  a variety  of  orbits  'he 
predictive  aspect  is  a very  important  one.  However,  the  data  used  to  make  a 
model  are  generally  several  years  old  before  thev  are  available  for  modeling  the 
environment.  Consequently,  by  necessity,  predictions  of  5 to  15  \ears  are  implied 
in  the  use  of  such  models.  This  predictive  aspect,  coupled  with  complexity  of  the 
natural  time  behavior  observed  in  the  particle  population,  place  a somewhat  unknown 
limit  on  the  accuracy  of  the  models  that  cannot  be  overcome  even  with  an  extensive 
measurement  program.  Although  detailed  historical  and  performance  data  exist 
for  companies  l.sted  on  the  New  York  Stock  Exchange,  the  prediction  of  the  price 
of  the  stock  a decade,  or  even  a year,  in  the  future  is  known  to  be  unreliable. 

With  this  analogy  in  mind,  we  will  now  outline  the  various  time  variations  that 
have  been  handled  in  the  production  of  the  models  to  date,  and  point  out  some  of  the 
difficulties  and  the  successes. 

Serious  modeling  of  the  radiation  belts  really  began  following  the  Starfish 
detonation  of  9 July  1962,  because  this  man-made  perturbation  had  serious  effects 
on  the  operation  of  a number  of  spacecraft  in  orbit.  The  inner  zone  electron  fluxes 
were  dominated  by  the  Starfish  residue  for  several  years,  which  is  illustrated  bv 
a comparison  of  the  radial  profiles  at  the  equator  of  AK-2  and  AK-5  ns  shown  in 
Figures  1-6  (pp  23-281. The  decay  of  these  electrons  was  observed  for  at  least 
8 years  and  it  was  necessary  to  model  this  decay  as  a fiinctton  of  L and  E in  order 
to  produce  AK-5  1975  projected,  where  observations  taken  during  1964-1966  were 
extrapolated  to  1975,  and  AK-6,  which  is  extrapolated  to  1980.  Typical  values  of 
decay  time  are  given  in  Figure  9 (p  44)*’  and  vary  from  50  to  370  days.  Effects  of 
the  Starfish  residue  had  to  be  taken  into  account  up  until  June  1970  in  small  regions 
of  the  inner  zone.  This  can  be  seen  in  by  referring  to  Figure  24  (p  55)*’  where  the 
cutoff  contours  for  the  Starfish  residue  are  shown  ns  a function  of  1.  and  K.  Con- 
sequently, both  current  models,  AK-5  1975  projected  (for  solar  minimum  condi- 
tions) and  AK-6  (for  solar  maximum  conditions),  for  the  inner  zone  have  been  de- 
rived bv  subtracting  a decaying  Starfish  residue  from  data  taken  more  than  10  years 
ago.  It  is  apparent  from  the  above  terminology  that  a long-term  solar  cycle  change 
for  the  inner  belt  fluxes  can  be  inferred,  at  least  for  electron  energies  below 
690  keV;  values  of  flux  change  ns  a function  of  K and  I are  given  in  Figures  86-93 


5.  Hilberg,  R.  II.,  Teague,  M.J.  and  Vetter,  J.l.  (1974)  Comparison  of  the 

Trapped  Electron  Models  AF-4  and  AF-5  with  AF-2  and  AF -3.  TJtTSHr 
74-13. 

6.  Teague,  M.J.,  and  Stassinopoulos,  F.  G.  ( 1972)  A Model  of  the  Starfish  Flux 

in  the  Inner  Radiation  Zone,  GSFC-X-601-72-4ft?^ 


(pp  134-141).  * The  flux  changes  are  typically  a factor  of  3;  however,  at  higher 
energies  no  such  change  can  be  obtained  because  of  the  masking  effects  of  mag- 
netic storms. 

The  effects  of  magnetic  storms  below  L 2.4  can  be  seen  in  Figures  8-11 
(pp  60-63).®  The  spectrum  hardens  during  the  storm  enhancements,  the  total 
energetic  electron  content  increases  by  factors  of  3 to  6 for  higher  L,  values  of 
this  region,  with  much  smaller  effects  in  the  heart  of  the  inner  zone.  These  en- 
hancements are  observable  for  2-3  months  and  are  definitely  discrete  events. 
Because  the  frequency  of  these  inner  zone  events  is  no  more  than  one  to  three  a 
year,  their  effects  cannot  be  handled  well  on  a statistical  basis.  The  qule;-time 
fluxes  with  solar  cycle  changes  have  been  used  along  with  an  average  contribution 
from  three  magnetic  storms,  so  that  the  average  inner  zone  flux  represents  a 
time-averaged  one  that  might  be  realized  over  a period  of  1 year.  If  one  were 
able  to  measure  the  omnidirectional  flux  at  one  point  in  space  and  integrate  this 
over  time,  the  fluence  would  be: 

T 

Ft(E,  D Do.  1.)  • J J(F,  B Bo.  1.,  t)  dt 
0 

and  the  time-averaged  flux  would  be: 

T - Kt/t  . 

The  omnidirectional  flux  of  AE-5  1975  projected  and  AE-6  corresponds  to  that 
which  would  be  expected  to  be  obtained  if  T were  approximately  1 year. 

Time  variations  for  outer  zone  electrons  require  a somewhat  different  treat- 
ment and  the  excursions  during  magnetic-substorm  enhancements  are  one  or  two 

orders  of  magnitude  larger  than  in  the  inner  zone  when  considered  as  a relative 

o 

change  in  the  ambient  flux.  A good  illustration  can  be  seen  in  Figure  9 (pp  41)' 
where  the  time  behavior  is  more  complex  at  l.  ■ 5.0  than  it  is  at  1.  4.0.  The 

approach  that  has  been  taken  in  the  outer  zone  has  been  to  use  a statistical  ap- 
proach to  describe  the  time  behavior  and  to  attempt  to  obtain  a time-averaged  flux. 
Here  the  averaging  time  needs  to  be  about  six  months  in  order  to  expect  the  observed 


7.  Teague.  M.J..  and  Vette,  J.l,  (1972)  The  Inner  Zone  Electron  Model  AF-5, 
NSSDC  72-10. 

8 Teague,  M.J.,  and  Vette.  J.l,  (1971)  Variation  of  the  Electron  Spectrum  in 
the  Inner  Radiation  Belt.  September  1964  to  Present,  \ssfk 1 71-1  1. 

9.  Singlev,  G.  \V. , and  Vette,  J.l.  (1972)  A Model  Environment  for  Outer  /.one 
Flections.  NSSDC  72-13. 
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fluence  to  approximate  that  computed  with  the  model.  Of  course,  the  local  time 
\ariation,  which  has  beer  discussed  earlier  as  a spatial  variation  must  be  removed. 
This  has  been  done  by  performing  a local  time  average  of  the  flux  at  a given  L 
value.  To  place  these  various  time  variations  into  perspective,  let  us  review  the 
situation  in  the  equatorial  geostationary  orbit  near  the  geomagnetic  equator 
(I,  » 6.6,  R B * 1).  The  instantaneous  excursions  can  be  up  to  three  orders  of 

O T 

magnitude  as  shown  in  Figure  11  (p  56)“:  this  data  set  is  fixed  in  local  time  at 
midnight  so  that  local  time  considerations  do  not  have  to  be  unfolded.  However, 
for  other  than  commensurntei y orbiting  spacecraft,  it  is  impossible  to  avoid  this 
unfolding  process.  As  we  have  remarked  in  Section  3,  the  local  time  variations  at 
L « 6.  6,  B B0  1 1 range  up  to  a factor  of  3.  2 at  1.9  MeV  and  larger  variations  can 
be  expected  at  higher  energies.  A graphical  example  of  this  can  be  seen  in  Fig- 
ure 22  (p  GO).  “ 

The  solar  cycle  effects  in  the  outer  zone  electrons  seem  to  be  confined  to  L 

values  below  5.0  and  must  be  inferred  by  studying  the  fluxes  averaged  over  about 

o 

one  year.  This  effect  is  illustrated  in  Figures  30  and  31  (pp  G2-63F  and  can  be 
seen  to  amount  to  maximum  changes  between  a factor  of  3 to  10,  depending  on  the 
energy  of  the  particle  and  the  L value.  Because  of  the  large  short-term  time 
variations,  it  has  not  been  possible  to  obtain  an  orderly  change  in  the  average 
fluxes  as  was  possible  in  the  inner  zone.  At  any  rate,  the  effects  of  solar  cycle 
changes  have  been  handled  by  producing  models  of  both  the  inner  and  outer  zone  at 
epochs  corresponding  to  minimum  and  maximum  conditions  (the  maximum  desig- 
nation is  only  approximate,  because  data  in  1969,  the  latest  solar  maximum,  were 
not  available  for  analysis  when  AE-4  and  AE-5  were  produced).  Considering  the 
magnitude  of  the  flux  changes  and  noticing  the  effects  of  magnetic  storm  enhance- 
ments in  the  inner  zone  and  the  local  time  effects  in  the  outer  zone,  a more  sys- 
tematic model  does  not  appear  warranted. 

Some  recent  work  by  Paulikas  and  Blake  (private  communication),  utilizing 
their  ATS  6 data  and  the  energetic  electron  data  of  Mcllwain  from  ATS  5,  demon- 
strate quite  clearly  that  at  energies  below  1.  55  MeV  there  are  no  solar  cycle 
effects  in  the  equatorial  geostationary  orbit  region.  This  is  illustrated  in  Figure  9. 
During  mid-1974  when  ATS  6 was  at  94°W  and  ATS  5 was  at  105°\V,  the  two  differ- 
ent experiments  were  intercompared  so  that  differences  of  energy  thresholds  and 
errors  in  geometric  factors  would  not  affect  the  results.  It  is  readily  seen  that  no 
systematic  solar  cycle  effect  can  be  ascertained  even  when  yearly  averages  are 
employed.  This  is  consistent  with  the  results  cited  earlier  where  the  same  exper- 
iment on  OGO  1 and  OGO  3 were  also  intercompared  in  the  same  time  period  by 
Winckler.  There  is  another  characteristic  that  can  be  noted  in  Figure  9;  the  yearly 
averages  for  the  electrons  above  3.  9 MeV  exhibit  a stochastic  change  that  is  about 
a factor  of  3.  Consequently,  an  inherent  error  in  these  more  energetic  electrons, 
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Figure  9.  ATS  5 and  ATS  fi  Daily  Flux  Averages  With  Sunspot  Numbers. 
(After  Paulikas  and  IMake,  private  communication) 
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even  for  yearly  averages,  must  be  anticipated  by  those  using  models.  New  datn 
will  be  required  to  determine  if  this  effect  is  larger  at  other  L values  and  at  H/13^ 
values  greater  than  1. 

Protons  exhibit  much  smaller  short-term  variations  than  do  electrons.  This 
is  amply  demonstrated  by  looking  at  Figures  1-6  (pp  23-26).  3 Certainly  the  more 
energetic  protons,  say  above  8 MeV,  show  very  small  changes  while  those  below 
1 MeV  can  show  changes  approaching  a factor  of  10.  There  have  been  two  obser- 
vations of  large  changes  of  protons  above  35  MeV  caused  by  magnetic  storms  that 
show  both  adiabatic  and  nonndinbatic  changes.  These  are  illustrated  in  Figures  1 
and  2 (p  11)*°  and  resulted  in  a reduction  of  the  particle  fluxes  which  recovered 
with  time.  The  bump  in  the  radial  profile  of  40  MeV  protons  discovered  by 
Mcllwaln  on  Kxplorer  15  nnd  Injun  3 and  studied  further  by  him  on  Kxplorer  26, 
provided  some  evidence  of  enhancements  of  energetic  protons  in  a confined  region 
of  space.  This  particular  distribution  discussed  by  Mcllwain  can  be  seen  In  Fig- 
ures 3 nnd  4 (p  11).  ^ Since  most  of  these  tyt>es  of  changes  are  not  predictable, 
occur  infrequently,  nnd,  for  nonndtabntic  chnnges,  require  a yenr  or  more  to  re- 
cover, they  are  best  handled  by  the  user  being  aware  that  such  departures  from  n 
stable  model  can  be  expected  on  rare  occasions. 

The  long-term  solar  cycle  behavior  at  low  nltitudes  has  been  taken  Into 
account  In  APB.  The  limited  region  of  H^H(1  I,  space  where  this  effect  is  observed 
can  be  seen  In  Figures  50-64  (pp  81-86),’’  These  changes  do  result  In  flux  level 
differences  of  more  than  an  order  of  magnitude,  but  the  altitude  gradients  are  so 
steep  in  these  regions  that  slight  differences  in  the  achieved  orbits  for  such  mis- 
sions could  result  in  equally  large  departures  from  the  prclnunch  expectations. 

The  time  variations  for  trapped  protons  do  not  seem  to  place  stringent  limits 
on  the  accuracy  of  the  fluxes  in  future  models.  However,  to  Improve  existing 
proton  models,  properly  calibrated  detectors  providing  complete  enough  coverage 
to  obtain  good  time-averages  will  be  required  to  reduce  errors  in  future  models. 


6.  Till*.  NKAR-TKRM  OUTI.OOK  FOR  KNI.ROI.TIC  PARTK:i.K  1'NVIRONMF.NTS 


The  electron  situation  in  the  equatorial  stationary  orbit  with  energy  range  up 
to  3.0  MeV  will  be  completely  defined  through  the  work  of  Pautikas,  Hlake,  nnd 
Hilton  (private  communication)  who  will  publish  a detailed  analysis  of  their  ATS  1 
nnd  6 results  along  with  intercompared  data  from  Mcllwain's  ATS  5 energetic 


10.  Lavinc,  J.P.,  and  Vette,  J.t.  (1070)  Models  of  the  Trapped  Itndintion 

Knvlronment- Volume  V:  Inner  licit  Prolons.  NASA  SP-3o!m. 

11.  Irvine,  J.P.,  nnd  Vette,  J.I.  (1070)  Models  of  the  Trapped  Itndintion 

Environment -Vohmv?  VI;  High  Hncrgy  ProTons,  SlASA  SP-3024. 


particle  experiment.  This  data,  along  with  that  from  OOO  5,  OV1-19,  A£l  R,  and 
ISIS  2,  is  in  the  process  of  being  analyzed  to  obtain  a new  outei  zone  model  AE-7 
This  effort  should  be  completed  in  about  a year.  Because  of  the  length  of  time 
necessary  to  produce  a properly  analyzed  model,  it  was  necessary  to  issue  Alll-7 
on  an  Interim  basis  to  account  for  the  new  electron  measurements  about  2 MeV, 

Some  of  the  changes  from  AM-4  that  can  be  expected  arc  illustrated  In  Figures  10 
and  11  where  the  comparison  of  AMI-7  with  some  data  are  shown.  The  HI  model  is 
based  on  the  OV1-19  analysis  of  Vnmpola.  As  one  can  see  in  Figure  10,  there  is 
basic  disagreement  between  the  AZl'R  and  OVl-19  results.  Both  measurements 
were  taken  in  the  same  time*  period. 

The  OV1-19  data  of  Vnmpola  is  presently  being  studied  along  with  West's 
OCiO  f>  data  to  resolve  the  discrepancies  with  the  AZUH  data  of  llovestndt.  Clearly, 
the  previous  All-4  model  is  low  at  energies  above  1. 5 MeV  based  on  the  ATS  li 
data  of  Paulikns  and  Blake  at  the  geostationary  orbit.  Since  ttie  calibration  of 
these  detectors,  as  shown  in  Figure  7,  have  been  carried  out  very  carefully  to 
high  energies,  the  validity  of  their  results  are  unquestioned.  In  addition,  their 

treatment  of  threshold  detectors  is  identical  to  that  which  we  introduced  in  AF-1 

2 

and  their  statistical  analysis  is  also  identical  to  that  which  we  introduced  in  the 
All-3  construction  and  used  in  making  All-4  (a  description  of  this  technique  can 
also  be  found  in  Ref  9).  The  major  change  from  All-4  is  a hardening  of  the  spectrum 
above  an  energy  dependent  on  I..  Whether  the  III  and  1X1  limits  used  in  AKl-7  will 
be  adopted  for  AM-7  is  not  known.  However,  with  the  stochastic  changes  observed 
even  in  the  yearly  averaged  fluxes  shown  in  Figure  0 at  3,9  MeV  such  a feature 
seems  justified.  Part  of  the  discrepancy  between  All-4  and  OV1-19  data  at  1 4.0 

may  In*  a solar  cycle  effect.  Solar  maximum  occurred  in  1009  and  the  latest  data 
used  in  AM-4  was  early  1907, 

Looking  beyond  the  present  analysis,  the  energetic  particle  data  that  will  be- 
come available  for  modeling  are  shown  in  Table  2 where  only  currently  operating 
spacecraft  are  given.  We  have  only  included  those  experiments  where  electrons 
of  at  least  0.8  MeV  and  protons  of  at  least  8 MeV  are  measured.  There  are  a 
number  of  other  experiments  on  these  spacecraft  and  MSA  GMOS  that  measure  lower 
energy  particles. 
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Table  2.  Operating  Spacecraft  With  Energetic  Trapped  Particle  Experiments 


I 


7 KK.COMMKMIATIONS  M>H  UTl  K!  Ml  AS!  HKMI  NTS 


In  order  to  produce  more  accurate  energetic  particle  environments  in  the 
future,  the  fallowing  recommendations  result  from  the  previous  discussion: 

(1)  Use  well  calibrated  omnidirectional  instruments  in  orbits  of  opportunity 
and  determine  saturation  effects  before  launch. 

(2)  Make  special  efforts  to  extend  the  electron  energy  range  up  to  at  least 
8 MeV:  in  doing  this  attention  should  In*  given  to  discrimination  against  the  con- 
tamination of  penetrating  protons  in  regions  of  the  inner  /one  where  their  fluxes 
might  be  large  relative  to  those  of  energetic  electrons. 

CD  Make  special  efforts  to  measure  protons  in  the  energy  range  from  100- 
t>00  MeV:  contamination  from  energetic  electrons  and  high  /.  particles  should  be 
discriminated  agninst  in  designing  such  systems. 

(4)  Use  digital  counters  and  a digital  spacecraft  data  system. 

(ft!  If  a single  orbit  for  trapped  flux  measurements  could  be  chosen,  select 
an  equatorial  orbit  with  about  1000  km  perigee  and  40,  000  km  npogce  altitude 
and  use  directional  instruments  that  can  obtain  the  complete  equatorial  pitch  angle 
distribution  over  the  energetic  pnrticle  spectrum.  A spin  ling  spacecraft  with  the 
spin  axis  normal  to  the  orbit  plane  is  ideal  in  this  case:  in  addition,  selecting  a 
12  hr  commensurate  orbit  might  greatly  aid  in  the  data  acquisition  coverage. 

(fd  Insure  that  resources  are  available  to  process  all  the  acquired  data  to  the 
point  where  engineering  corrections  are  made,  saturation  effects  are  removed, 
periods  of  suspect  instrument  performances  nnd  bad  data  are  deleted,  and  orbit 
(also  attitude  for  directional  instruments)  data  are  merged  with  the  accumulated 
corrected  counts. 

(7)  Insure  that  proper  documentation  nnd  calibration  curves  are  provided  to 
the  environment  modeling  group  along  with  the  reduced  data  within  l year  nfter 
launch.  This  will  permit  more  timely  nnd  accurate  models  to  lie  produced. 

(8)  Insure  that  orbital  position  accuracy  is  adequate  to  handle  measurements 
pro|ierly  where  the  spatial  gradients  are  large.  These  regions  are  at  low  altitude 
and  in  the  I.  range  between  2 nnd  4 where  protons  exhibit  steep  radial  and  energy 
gradients. 

(0)  If  H/ll  nnd  I coordinates  are  computed,  insure  that  the  same  model  of 
the  internal  field  is  used  nnd  appropriate  secular  changes  are  incorporated  so  thnt 
comparisons  of  measurements  from  different  satellites  can  be  conducted  with 
minimum  error. 
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